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\  t  nn  .'  i.'-ioual  molocular-dvnamics  model  of  a  chemically-sustained  shock  wave  in  a  molecular  solid  is  dis 
(  us'e.i  A  nody  expression  that  realistically  describes  chemical  bonding  in  condensed  phases  is  used  for  the 

interatomic  potential.  1  he  results  demonstrate  that  simple  reactive  collisions  are  sufficient  to  self-sustain  a  shock 
wave  with  a  velocity  that  is  comparable  to  experimental  detonation  velocities.  Features  characteristic  of  macroscopic 


detonations  such  as  an  intrinsic  detonation  velocity  and  a 
l.  INTRODUCTION 

Molecular  dynamics  (MD)  simulation  has  proven  to 
be  an  excellent  method  for  studying  atomic-scale  reac¬ 
tion  dynamics'  in  a  of  variety  systems.  For  this  technique 
to  be  of  use,  however,  the  system  to  be  studied  should 
be  restricted  to  ato  nic  dimensions  and  short  timescales 
(up  to  nanoseconds),  and  an  adequate  potential  energy 
surface  should  be  used.  Nonreactive  shocks  provide  ideal 
phenomena  to  be  studied  using  MD  techniques;  the  re¬ 
sults  of  MD  studies  using  Lennard- Jones  potentials  have 
compared  very  well  to  macroscopic  models,1  suggesting 
that  shock  waves  in  liquids  and  solids  can  be  well  de¬ 
scribed  on  an  atomic  scale.  Detonating  solids,  however, 
require  some  form  of  energy  release  to  self-sustain  the 
shock  wave,  and  so  Lcnnard-Jones  potentials  are  gen¬ 
erally  inadequate.  Furthermore,  to  realistically  model 
chemical  energy  release  the  forces  should  be  many-body 
in  nature  since  the  reactivity  of  a  given  atom  depends  on 
its  bonding  to  surrounding  atoms.  Analytic  expressions 
that  realistically  describe  changes  in  chemical  bonding  in 
the  solid  state  have  only  recently  been  introduced,  and 
with  few  exceptions2  have  been  restricted  to  metals  and 
semiconductors.3 

In  this  work  we  present  a  two-dimensional  micro¬ 
scopic  model  for  a  detonating  molecular  solid  that  is 
based  on  the  many-body  Tersoff  covalent  bonding  for¬ 
malism.*  The  simulations  demonstrate  that  macroscopic 
features  of  detonating  solids  such  as  an  intrinsic  detona¬ 
tion  velocity  and  a  following  flow  are  apparent  on  the 
microscopic  level.  The  simulations  display  a  detonation 
front  that  is  on  the  order  of  a  few  atomic  dimensions 


following  flow  are  observed  on  a  microscopic  scale. 

wide  and  a  high  density  reaction  region  where  much  of 
the  molecular  nature  of  the  solid  is  lost. 

II,  MODEL 

TersofT  has  recently  introduced  a  simple  analytic  ex¬ 
pression  that  realistically  describes  bonding  in  covalently- 
bonded  solids.*  Encouraged  by  this  work  (and  additional 
extensions  of  this  approach  to  few-body  reactions5),  we 
have  adapted  a  simplified  form  of  the  TersofF  expression 
to  model  an  energetic  molecular  solid.  In  this  approach 
the  binding  energy  is  written  in  a  form  similar  to  a  pair 
potential: 

&  =  £  £  ■  (U«(r'r)  -  Bij  '  VB(r„) 

•  ;(>•) 

-O-0.j)-V„B(r.j)l.  (1) 

The  functions  Vg(r ),  Vg(r)  and  Vfjg(r)  represent  repul¬ 
sive,  attractive  bonding,  and  nonbonding  pair-additive 
interactions,  respectively,  r,;  is  the  scalar  distance  be¬ 
tween  atoms  i  and  j,  and  /c(r)  is  a  function  that  re¬ 
stricts  the  range  of  the  potential.  The  many-body  aspect 
of  the  potential  enters  through  Bij  which  weakens  the 
pair-additive  bonding  term  (and  hence  the  bond  energy) 
between  atoms  t  and  j  as  the  local  coordination  of  the 
atoms  increases.  For  our  model  it  is  also  used  to  modu¬ 
late  the  nonb^nded  interactions. 

Following  Tersoff,  the  function  BtJ  is  given  by 

Bn  =  (Bn  +  Bn)/ 2  (2a) 


/A,  (I  :  <•  1  l>\"l  (r,f  -  1  ,iL  )]  }  (2/*) 
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(  j)  III  RESULTS  A\n  DISCI  SSION 


arc  used  for  tin*  repulsive  and  attractive  bonding  pair 
interactions,  the  uonbonded  pair  term  is 

l  -  (■->) 


and  fc(r)  is 
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A  Di  <  r  <  D2 

r  >  Di. 
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Eqs.  1-6,  while  maintaining  the  simplicity  of  a  pair 
potential,  incorporate  general  bonding  concepts.  In  the 
limit  of  isolated  diatomic  molecules  each  Bi}  equals  one 
and  F,q.  1  reduces  to  a  sum  of  morse  functions.  As 
atomic  coordination  increases,  the  pair  additive  bond 
strengths  decrease  and  the  equilibrium  distances  increase 
in  a  way  that  yields  the  Pauling  bond  order  relationship.4 
Since  the  functional  forms  are  relatively  simple,  this  ap¬ 
proach  is  efficient  for  modeling  the  simultaneous  reaction 
of  many  molecules,  such  as  occurs  in  a  condensed-phase 
detonation. 


Rather  than  try  to  describe  one  particular  system,  we 
have  developed  a  simple  generic  model  for  an  energetic 
molecular  solid.  The  parameters  entering  the  function 
B,,  are  chosen  so  that  bonding  pair  terms  are  signifi¬ 
cantly  weakened  for  atomic  coordinations  greater  than 
one.  This  insures  that  the  lo'vcst-encrgy  state  of  the  sys¬ 
tem  is  a  collection  of  diatomic  molecules,  and  maintains 
the  molecular  nature  of  the  solid  at  thermal  tempera¬ 
tures  and  low  pressures.  They  have  also  been  chosen  to 
reproduce  chemically-sound  potential  barriers  for  reac¬ 
tive  collisions.  Energy  release  is  incorporated  by  using 
,  two  different  kinds  of  atoms  with  different  well  depths  in 
;  th/c  morse  functions  (Eqs.  3  and  4).  For  hetcronuclear  di¬ 
atomic  molecules,  a  well  of  depth  of  2  eV  is  used  while  for 
homonuclear  diatomic  molecules  the  well  depth  is  5  eV. 
The  initial  unreacted  system  is  composed  of  heteronu- 
clcar  diatomic  molecules  that  are  bound  into  a  solid  by 


The  shocks  were  initiated  by  striking  a  two  dinien 
siona!  crystal  with  a  simulated  flyer  plate  that  was  com 
posed  of  the  same  material  as  the  initial  crystal,  l’oriodh 
Loju.laries  were  maintained  in  i  he  diieclion  peipmdu  u 
lar  to  the  shock  Shown  in  Figs,  la  and  lb  ate  snapshots 
from  a  simulation  where  the  flyer  plate  was  given  an  ini¬ 
tial  velocity  of  6  kui/dct.  in  i  ig.  la  tin.  shock  car, 

be  identified  by  a  distortion  of  the  crystal  structure.  A: 
about  10  A  behind  the  plate-induced  shock  front  some  of 
the  molecules  have  collided  and  begun  reacting.  Fig.  1!> 
shows  the  system  after  «S  ps.  At  this  point  the  initial 
shock  caused  by  the  flyer  plate  has  diminished,  while  the 
chemically-sustained  shock  continues  to  propagate.  A 
hot  gas  consisting  mostly  of  product  molecules  and  a  few- 
free  atoms  and  reactant  molecules  is  visible  expanding 
out  the  back  of  the  shock  at  the  left  of  Figs,  la  and  lb. 
and  a  high  density  region  separating  the  gas  from  the 
shock  front  is  visible  in  the  center  of  Fig.  lb. 

The  velocities  of  the  front  edge  of  chemically-sustain¬ 
ed  shocks  versus  time  are  given  in  Fig.  2  for  two  different 
simulations.  The  dotted  line  is  from  the  simulation  de¬ 
scribed  above  and  the  solid  line  is  from  a  simulation  that 
used  an  impact  plate  velocity  of  10  km/scc.  In  the  latter 
case  the  chemical  reactions  began  at  the  plate-solid  in¬ 
terface,  with  no  leading  platc-induced  shock  wave.  After 
a  brief  period  both  of  the  chcmically-drivcn  shocks  reach 
the  same  constant  velocity  of  4.8  km/sec.  The  steady- 
state  velocity  of  the  chemically-sustained  shock  is  appar¬ 
ently  an  intrinsic  property  of  the  microscopic  system. 

Despite  the  microscopic  size  of  the  system  simulated 
here,  the  model  displays  behavior  characteristic  of  macro¬ 
scopic  detonations.  First,  for  the  5  km/sec  impact  plate 
velocity  a  nonreactive  compressional  wave  preceeds  the 
detonation  front.  This  is  similar  to  real  systems  where 
detonation  is  initiated  behind  a  leading  shock  wave.0 
Second,  the  detonation  velocity  predicted  by  our  sim¬ 
ple  model  agrees  very  well  with  experimental  velocities 
(typically  3-7  km/sec).  It  should  be  emphasized  that 
no  adjustment  of  the  potential  was  made  to  achieve  this 
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FIG  1.  Snapshots  from  the  simulation.  The  two  types  of  atoms  are  drawn  with  different  sizes,  and  lines  connect 
atoms  that  are  within  1.4  A.  The  shock  is  traveling  from  left  to  right,  a)  The  system  from  -20  A  to  100  A  at  1  ps. 
b)  The  system  from  2S0  A  to  300  A  at  S  ps. 


velocity,  rather  it  is  a  prediction  of  the  model.  Further¬ 
more,  this  velocity  is  intrinsic  to  the  material,  again  in 
agreement  with  macroscopic  systems.  Third,  the  high- 
density  region  (apparent  in  the  center  of  Fig.  lb)  exhibits 
a  steady  forward  flow  of  «4.2  km/sec,  and  its  width  in¬ 
creases  linearly.  For  long  times  this  region  would  presum¬ 
ably  evolve  into  a  macroscopic  following  flow,  in  agree¬ 
ment  with  continuum  models.7 


Time  (ps) 

FIG.  2  Velocities  for  chemically-r  istained  shock  fronts 
versus  time. 


The  atomic-scale  structure  and  dynamics  of  the  chem¬ 
ically-sustained  shock  is  of  particular  interest,  since  this 
regime  in  a  real  detonation  is  not  easily  accessible  to  ex¬ 
periment.  Shown  in  Fig.  3  is  the  potential  energy  per 
particle  versus  position  (relative  to  the  shock  front)  cal¬ 
culated  from  the  simulation.  The  unreacted  material  is 
to  the  right  of  the  figure  and  the  shock  front  is  at  zero. 


Distance  (A) 

Fig.  3.  Potential  energy  versus  position  relative  to  the 
front.  The  front  is  at  0  A  and  unreactcd  material  is  to 
the  right.  Expansion  into  the  gas  phase  begins  at  ssSO  A 
behind  the  front. 
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■  *  p-  >:,!<■  I  In-  Miuiil.il  urn  At  I  lie  shock  iioiil  ilir 
potential  rneigy  uses  sharply  to  a  value  that  remains 
d.tislaiii  ihoughout  tin*  high  density  legion.  att<l  the  po¬ 
tent  ial  <  :  does  nut  drop  to  a  value  lower  I  itait  t  lie  till 

ie.i'  led  mateiial  until  the  molecules  expand  into  the  »,is 
phase.  Apparently  chemical  energy  is  not  fully  accessed 
until  the  gas  phase  molecules  are  formed.  This  picture 
of  the  potential  energy  corresponds  well  with  Fig.  )h: 
although  dirt  met  molecular  units  are  apparent,  as  reac¬ 
tant  and  product  molecules,  in  the  high-density  region 
the  molecules  are  undergoing  multiple  collisions  and  the 
molecular  nature  of  the  solid  is  lost. 

Several  extensions  to  the  model  are  possible.  First, 
only  one  choice  of  parameters  for  the  potential  energy 
expression  has  been  discussed.  The  functional  form  is 
sufficiently  flexible  that  features  of  the  potential  energy 
surface  can  be  changed  so  that  relationships  between  re¬ 
action  kinetics  and  detonation8  can  be  studied.  Second, 
the  microscopic  characteristics  of  the  simulation  may  be 
influenced  by  the  reduced  dimensionality  of  the  simu¬ 
lations.  Because  the  form  of  the  potential  is  not  re¬ 
stricted  to  low  dimensions,  extension  of  our  model  to 
three-dimensions  is  straight  forward.  Third,  although 
changes  in  bonding  are  realistically  described  by  the  mod¬ 
el,  the  system  presented  uere  is  simple  compared  to  real 
energetic  materials.  More  complicated  forms  of  the  po 
tential  energy  expression  given  above  have  been  used  to 
decribe  condensed- phase  carbon9  and  the  exchange  reac¬ 
tions  H+H2  and  0+02,5  and  so  the  further  development 
of  our  model  to  describe  real  systems  appears  promis¬ 
ing.  Finally,  no  attempt  has  been  made  to  include  ef¬ 
fects  outside  of  simple  nuclear  collisions.  Although  such 
effects  may  ultimately  prove  important,  the  results  pre¬ 
sented  here  suggest  that  simple  collisional  models  are  a 
good  starting  point  for  describing  many  of  the  properties 
of  detonating  solids. 
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